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Abstract 

We propose an algorithm which produces a randomized strategy reaching optimal 
data propagation in wireless sensor networks (WSN). In jHj and [S], an energy bal- 
anced solution is sought using an approximation algorithm. Our algorithm improves 
by (a) when an energy-balanced solution does not exist, it still finds an optimal so- 
lution (whereas previous algorithms did not consider this case and provide no useful 
solution) (b) instead of being an approximation algorithm, it finds the exact solution 
in one pass. We also provide a rigorous proof of the optimality of our solution. 

1 Introduction 

In |>] , the problem of finding and energy-balanced solution to data propagation in a Wire- 
less Sensor Network (WSN) is considered. They model a WSN by splitting in it in slices, an 
try to maximize lifespan of the network by ensuring that the expected energy consumption 
in each slice is the same. The case they consider is the following: sensors are distributed 
randomly and informally in a circle (or a sector), and data has to be propagated from the 
WSN towards the sink, which is in the center of the circle. It is considered that data is 
to be slided towards the sink along the network in a multi-hop fashion, but that any node 
which contributes to this propagation can decide to eject the message directly to the sink 
in an energy expensive single-hop long range data transmission. Using this model, the 
authors of jH] show that energy-balance can be reached if a recurrence relation between 
the probabilities that a slice ejects a message single-hop to the sink is satisfied. They then 
propose an approximation algorithm (converging exponentially fast towards the solution) 
to compute the solution to the energy-balancing problem. Although they do not prove it, 
it happens to be that this energy-solution always exists in the context they have chosen: 
a circle (or a sector) with informally distributed sensor nodes inside of it, and a uniform 
distribution of events (which generate data to be slided towards the sink) inside the WSN. 

*{pierre.leone|olivier.powell|jose.rolim}@cui.unige.ch 



1 



In 0, a more general case is studied: events may not happen according to a uniform dis- 
tribution, and the sensors may not be distributed uniformly over the area to be monitored, 
Moreover it is not assumed that the distribution of events is known a priori. A solution 
to the energy-balance problem is then computed on-line (by the central sink) whilst the 
distribution is discovered from observation. They show that their algorithm tends to the 
energy-balanced solution (when this solution exists). However, neither [S| nor [H] consider 
cases where an energy-balanced solution does not exist. The idea we want to fulfill is 
to find an optimal solution, whatever the distribution of events and the distribution of 
sensors is. We provide an algorithm to compute such a solution, which turns out to be 
the energy-balanced solution when this solution exists. Another improvement is that our 
algorithm, in the contrary to the algorithm provides an exact optimal solution, whereas 
the algorithms from [B] and [S] are approximation algorithms converging exponentially, 
but with potential infinite run time. Our algorithm is based on a very simple idea, which 
can informally be stated as if your neighbour has more energy available then you, let him 
do the job. Although this is over-simplifying and special cases have to be treated for the 
general solution, we propose as an application of this idea a completely distributed solution 
(i.e. not computed in a centralized sink) which converges to the optimal- solution under 
some reasonable assumptions, and provide some simulation evidence of the interest of this 
approach. 

2 Model 

The model we study is the same as the one from jHj and [Hj, and is the following. A sensor 
network is composed of iV slices {Si} 1<i<N . Each of these slices represents a subset of the 
sensors of the network. For 2 < j < N — 1, the jth slice, Sj, has two neighbour slices: 

1. Sj+i, which is called the previous slice (with respect to Si). 

2. Sj-ii which is called the next slice. 

The first slice, Si, is also neighbour to a special node of the network, called the sink. The 
last slice, Sn has no previous slice, but only a next slice: Sjv-i- 

Definition 1 We say that a slice detects an event if one of the sensors of the slice detects 
an event. 

When a slice Sj detects an event, it needs to report by sending a message to the sink. 

Definition 2 A generated message is a message which has to be sent by a slice Si following 
the detection of an event by Si. 

When a slice is in presence of a generated message, it has the choice between two strategies: 

1 . Send the message directly to the sink: we say that the (generated) message is ejected 
by Si. 
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2. Send the message to the next slice, Si-i, putting the responsibility of communicating 
the message to the sink on we say that the (generated) message is slided along 
the network from Si towards Si-i. 

Definition 3 Once a generated message has been slided along the network, we do not call 
it a generated message anymore, but we rather call it a sliding message. 

When a slice receives a sliding message from its previous slice, it has the choice between 
the two same strategies as for generated messages: either send the message to its next slice 
(the message keeps sliding along the network towards the sink) or stop the message, and 
send it directly to the sink: we say that the sliding message is ejected to the sink. With 
this model in mind, we describe a slice and its behaviour by the following parameters. 

• gi is the number of events being detected by Si. 

• hi is the amount of energy available in slice Si, which is proportional to the number 
of sensors, assuming each sensor has the same amount of energy available. 

• di is the distance of Si to the sink. 1 

• fi is the number of messages slided from Si towards Si-±. 

• ji is the number of messages ejected from Si to the sink. 

• Pi is the sliding probability: Pi = j^p- 

• €i is the ejection probability by slice Si of a sliding message. 2 
Notice that the following recurrence relation holds 

fi + ji = fi+i + 9i (1) 

Transmitting messages has a price, which has to be payed by the transmitter. In this 
model, we take the following convention. 3 

Convention 4 Sliding a message from Si towards S^-i costs 1 energy unit to Si, and 
ejecting a message from Si to the sink costs d\ energy units to Si . 

lr rhe distance of a slice to the sink is the maximum over all distances between sensors of the slice and 
the sink. We suppose w.l.o.g. that d\ > 1 for every i, and that d i+ i > di. 

2 That is, ej is the fraction of the /j+i sliding messages coming from Si+\ which are ejected from 5, to 
the sink. E.g. if gi = 0, then a = 1 — pi. 

3 With this convention, we do not consider the cost of receiving a message. However, this does not yield 
much loss of generality. For example, if we wanted to fix the cost of receiving a message to be equal to 
short range transmission (i.e. 1 energy unit), we could include this in the model by setting the price of 
long range transmission to df + 1 and the price of short range transmission to 2. 
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Therefore, the mean energy consumption of slice Si is equal to the following: 

I'; ■ fi ■ J,<lj (2) 



Definition 5 We say that the network is energy balanced if for 2 < N it holds that 

1 1 

-r-Ei = - — Ei_ x (3) 

Oi Oi-i 

That is, the network is energy balanced if the mean energy consumption is the same for 
every sensor of the network. We next want to introduce a measure of the lifespan of a 
network pursuing the idea that the network is down as soon as any single slice runs out of 
battery. The lifespan of the network thus depends on the worst case energy consumption 
amongst slices. 



Definition 6 The lifespan of the network is defined as min {jf: j 



Ki<N 



Let m be a message either generated in or slided to the slice S^ Suppose Si ejects the 
message directly to the sink with probability q. The mean energy consumption per sensor 
of slice Si induced by the treatment of message m is equal to 

£[(1 -?)+?<$] ( 4 ) 

Since message m is slided towards Si-i with probability 1 — q, it generates an average 
energy consumption per sensor of Si-\ equal to the following: 

-L(l-g)[(l _£._,) + e l _ 1 d 2 l _ 1 ] (5) 

0;-i 

since is the ejection probability of a sliding message for S^-i. Therefore, we get the 
following claim. 

Claim 7 Let m be a message that has to be treated by Si. Suppose that m is slided (by Si) 
along the network towards Si-i with probability 1 — q and that for for all 2 < i < N the 
following recurrence relation holds: 

_L [(i _ e .) + tid J\ = (i _ e.) J_ [(i _ ei _ ± ) + e^d^] ( 6 ) 

Oi Oj-l 

Then if q = 6j the treatment of message m (which implies contribution by all slices Si, S^i, 
. . ., Si) generates the same mean energy consumption per sensor in all slices Si, . . . , S±. 
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This follows from the fact that treating m has an average cost per sensor of slice Si (when 
q = €i) of m = (1 — + e^ 2 ). The mean average cost per sensor of is then seen 
to be equal to m too, because of (jHJ) and (JHJ), and so on for Si-2 to S\ (details left to the 
reader) . 

Relation (jSJ can be rewritten in the following form: 

ggi (1 + 6^ - 1) - 1 

^ (1 + e^ 2 - 1) + d m 2 - 1 U 

Setting t\ — \ (which is natural since the first slice can only send messages directly to the 
sink), we can compute directly the e$'s satisfying the recurrence relation for fixed b^s and 
di's. We next show how to use this fact to find an optimal solution for the problem of 
lifespan maximization. 



3 The Algorithm 

The algorithm we propose is an algorithm which computes a randomized strategy for data 
propagation in a sensor network separated in N slices of nodes at increasing distances 
from the sink. The input is a description of the network and a statistical description of the 
data to be propagated in the form of three sequences of same length {bi} 1<i<N , {di} 1<i<N 
and {gi} 1<i<N where the 6,'s describe the energy available in each slice, the diss's are the 
distances of each slice to the sink, and the g^s are the distribution of events generating 
data to be propagated in the network (c.f. section EJ. The output is a sequence {pi} 1<i<N 
representing a randomized strategy the network should apply in order to maximize its 
lifespan: each slice Si for 1 < i < N, when treating (either a sliding or a generated) 
message should send the data directly to the sink with probability 1 — pi and slide it along 
the network to the next slice with probability Pi. We can now start the description of the 
algorithm 

Since the frj's and dj's are given from the input for 1 < i < N, we can compute the e^s 
( 1 < i < N) which are the solution to relation © with e% — 1. 

Remark 8 The attentive reader may have noticed that although it holds that equation 
(0) implies that 6j < 1, nothing guarantees that t{ > 0. For simplicity, let us make the 
temporary assumption that ti > 0. We treat the case with negative e's in section Y3.1. 6 A 

For each 1 < i < N, we define the following: 

1. Gi is the number of messages to be treated at slice Si. 

2. Fi is the number of messages forwarded from slice S^+i towards slice Si. 

3. Ji is the number of messages ejected from slice Si directly to the sink. 

4. Ei is the average energy spent by slice Si, which is equal to Fi + Jidj 2 . 
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Using the g^s from the input we initialize F{ = Ji = and Gi = gi for every 1 < % < N. 
We then treat each slice one at a time, in a top-down fashion (i.e. from Si, close to the 
sink, towards Sn, away from the sink ): 

• Let Ji — G\. That is, all the messages generated at 5*1 are ejected to the sink. This 
means an average energy consumption in slice Si of E\ = G\d\ = Jidf + F\. 

• To take account of the fact that all the Gi messages to be treated have been ejected, 
we update the value of Gi to G\ = 0. 

We then repeat the following for each of the slices Si for i from 2 to N: First, let Jj : = 
bi a Ei-i . This means that the following holds: 

T-Jidi 2 = - — Ei-x (8) 
h 6i_i 

The interpretation of equation (jHJ is that the average cost per sensor in slice Si due to 
ejection of messages is equal to the average energy consumption per sensor (so far) in slice 
Si_i. So at this points, the average energy consumption per sensor in slice Si is equal 
to the average consumption per sensor in slice Sj_i (and, by induction, in Sj_ 2 , Sj_ 3 , • • • 
). Notice that if Jj > Gi, it means we are trying to eject more messages than the total 
amount of messages available to be treated, which is not physically possible. Therefore, for 
our approach to work smoothly it should be that there are enough messages to be treated 
at slice Si, i.e. it should be that the following holds: 

Ji = - — Vl^-i — (9) 
bi-idi 

This is an important condition, which does not hold without loss of generality. We explain 
how to overcome this limitation in section 13.1.11 For the time being, we consider only 
the case where the initial Gj's are large enough so as to ensure that equation Q holds. 
Although so far energy balance is reached for the slices Si to Si (because of equation (jHJ) 
and by induction), there are still Gi — Ji messages generated at slice Sj which have to be 
treated. Let us update Gi to Gi := Gj — Jj, thus Gi now counts the number of messages 
generated at slice Gi yet to be treated. These messages have to be treated in such a 
way that they will increase the average energy spent per sensor in each of the slices Sj 
to Si by exactly the same amount. The strategy is the following: a fraction (1 — q) of 
the Gi messages yet to be treated are forwarded to the next slice Sj_i, while the rest is 
ejected directly to the sink, thus increasing the average energy spent per sensor at slice Sj. 
Formally, this means setting the following: 

Fi := (l-€i)Gi (10) 
Jj := J i + t l G l (11) 

Of this fraction Gi (1 — Cj) of messages slided from Sj towards the slice Sj_i, a fraction 
1 — Cj_i will be further slided from Sj_i towards Sj_2, while the rest is ejected directly to 
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the sink from Sj_i, etc... Algorithmically, this means doing the following: 



Ft. 
Ji- 
Fi- 

Ji- 

F- 



-1 




Fi-i + (1 — (1 — e*) Gi 


(12) 


-1 




Ji-i + Q-i (1 — e») 


(13) 


-2 




Fi_ 2 + (1 - Ci _ 2 ) (1 - ei-x) (1 - ei) d 


(14) 


-2 




Ji-2 + ei_2 (1 - (1 - e») ^ 


(15) 


-3 




Fi_ 3 + (1 - e^) (1 - ei- 2 ) (1 - Q-i) (1 - ei) ^ 


(16) 


-3 




+ ei-3 (1 - ei-2) (1 - (1 - e,) 


(17) 



'1? 



During this phase 4 (forwarding/ejecting from slice to slice), the energy consumption per 
sensor at slice Si has been increased by := (1/6$) G^ [(1 — e$) + ejdf], the increase of en- 
ergy spent per sensor at slice Sj_i is equal to mj_i := (l/6j_i) Gi (1 — e$) [(l — 6j_i + ei_idf_ 1 )] . 
etc... But because the e^s satisfy equation (jZJ), all these m/s (1 < j < i) have the same 
value, as follows from claim So when we finish treating slice Si the average energy con- 
sumption per sensor in Si, £2, ... , S* is equal (slightly increased from the previous step), 
etc... We then go on to treating the next slice (£j+i) until we reach the last slice, Sjy. At 
this stage, the average energy consumption per sensor in each of the slices will be equal, 
thus energy balance will be reached. At this point and for each 1 < i < N slice Sj treats 
a total of Fi + J t messages, of which Fj are being slided and Jj are being ejected. The 
output of the algorithm representing the randomized strategy the network should apply is 
therefore going to be the following ordered sequence: 



Fi+Ji 



Ki<N 



3.1 Special cases 

In this section, we lift the restrictions following from the assumption that all e's are positive 
(remark |HJ), and from the assumption that equation Q holds, starting with the former. 



3.1.1 First Case: Little Messages 

Suppose that while executing the algorithm from previous section, equation Q does not 
hold for some i. What this means is that for some i, while treating slice Si, even if all the 
gi generated messages are ejected to the sink, Si will not spent as much average energy per 
sensor as slice Sj_i does; informally: there are not enough generated messages at slice Si 
(thus the title of this subsection). So what can we do? In essence, the solution is to get 
slices previous to Si (i.e. Sj + i, Sj+2, etc..) to forward some of there generated messages 
towards Sj, so that S, can "catch-up" with Sj_i. Practically, this is done by recursively 
using the algorithm described erstwhile. To bring this idea to effect, we go down one-level 
in the recursion. This brings the need to stack some values, which we will be able to 

4 Also, after this we can set G, = 0, since all messages at slice Si have been treated. 
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unstack when coming up one level in the recursion, and to compute some new values. The 
following enumeration stepwise explains how to implement this scheme. 

1. We start by stacking previous values and computing new ones. 

(a) Stack the current value of a variable start, which remembers at what position 
the last recursion started. If this is the first level of recursion, we stack the value 
start = 1. The new value of start is set to start = i, which is the position of 
the slice which is trying to catch-up. 

(b) Stack the current value of a variable maxNrg. If this is the first level of re- 
cursion, we stack the value maxNrg = oo. The new value of maxNrg is set 
to 

maxNrg = (F star t-i + J ^start-id start— l) 

Ostart-1 

which is the amount of energy to be spent by S start in order to "catch-up" with 

S start-l- 
ie) Stack the value of the previous e's. Set new current values for e^-'s for start < 
j<N: 

• e star t is set to 1 (slice S sta rt ejects every sliding messages it sees, in order to 
try catching- up with S sta rt-i)- 

• The other e^'s with start < k < N are computed using equation (|7J). 

2. Next, we go down one level in the recursion, which essentially means redoing the 
algorithm from section El but using the above newly computed e's, and considering 
the possibility of either going further down one recursion level, or on the contrary 
coming back up one recursion level. Before describing in more details what our 
algorithm does whilst going up/down recursion levels (as described above), let us 
make the following remark. 

Remark 9 While in the lower level of recursion and while treating slices S sta rt+i, S star t+2, 
etc., three cases may occur (for the sake of comprehensiveness, suppose we are treat- 
ing slice Sk for some k > 0): 

• Slice S s tart may be able to catch-up (using messages previously slided from slices 
Sstart+i t° Sk-i and the newly messages slided from Sk), i.e. sufficiently mes- 
sages will have been slided from its previous slices. In this case, we shall have 
to go back up one recursion level. 

• Sk may in turn not be able to spend as much energy as Sk-i, in this case, we 
need to further go down another level. 

• Finally, neither of the two cases above may happen, slices Sk, Sk-i, ■ ■ ■ , S sta rt are 
spending the same amount of average energy per sensor (but less than S star t-i), 
and we start treating the next slice, Sk+\ 
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We now describe in more details the way the algorithm from section |3] should be 
adapted in order to implement the above remark. Suppose we are in the first recursion 
level. We are treating slice Sk for some k > start. Slice S star t is trying to catch-up 
with slice S sta rt-i, and all the slices from S sta rt to Sk-i are spending the same average 
energy per sensor (which we can suppose to be true by induction). 

3. // there are not enough messages for S& to spend as much energy as Sk-i, i.e. if 
-^Gkd\ < y— (Fk-i + Jk-id\_]) we just do our best: we eject Jk = Gk messages, 
set Gk = 0, and further go down one recursion level. Now S& will have to try to 
catch-up with Sk-i- 

4. Else, we eject J k = -Jy — (Fk-i + J/ C _i(i|_ 1 ) messages from Sk and set Gk = Gk — Jk- 

k fe— 1 

Sk now has the same mean energy consumption per sensor as Sk-i, but there are still 
Gk messages to be treated. 

Remark 10 Essentially, we now want to slide the remaining Gk messages along the 
network, from Sk to S s t ar t, but with some precaution: 

• First of all, we still have to take into the account the e's and eject "e fractions" 
of the messages sliding along the network from Sk to S s t ar t- This was explained 
in detail in equations il(J\) to \18\). and we therefore we consider from now on 
this to be implicit whenever we talk of sliding messages along the network. 

• If we have enough remaining messages, i.e. if Gk is still sufficiently large, we 
will be able to let slice S start (and slices Sk to S s t ar t+i) catch up with S s t a rt-x- If 
this is so, we want to slide just enough messages so as to catch-up, then go back 
up one recursion levef , and only then shall we slide the eventually remaining 
messages along the network. 

Here is how we propose to implement the above remark. 

5. Set A nrg = maxNrg - £ (F k + J k d 2 h ) and msgToGolIp = j£f ( " r _f ) • Let $ = 
max{Gfc, msgToGolIp}. From the remaining Gk messages, we further slide and 
eject respectively F — (1 — e&) • $ and J — 6% ■ $ messages 6 . Finally, we need to make 
the following adjustments: F k = F k + F, Jk = J start + J and G = G — F — J. 7 

6. If there are enough messages for S s tart to catch-up with S s tart-i we can go back up 
one recursion level: 

This is to ensure that we can unstack the appropriate e's while coming-up a recursion level. E.g. there 
is no more reason to use the e's where e s tart was forced to the value e = 1 if S s t ar t has managed to catch-up 

with S s tart-1- 

6 This ensures that slices Sk to S s tart still have the same mean energy consumption per Sensor and that 
this is bounded by the mean energy consumption per sensor of slice Sstart- 1 

7 The Fj's and Jj's for start < j < k also have to be adjusted, acknowledging the fact that F new 
messages are slided along the network from Sk, and using the e's, but as explained in the first point of 
remark ITUl we consider this to be implicit. 
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• Unstack the previous e's. 

• Unstack the previous value of maxNrg. 

• Unstack the previous value of start. 

7. • If there are no more messages to treat for slice Sk, i.e. if Gk = 0, we can start 
to treat the next slice, Sk+i- This means jumping to point El above, but this 
time with k = k + 1. 

• Else we need to treat the remaining Gk messages. This is done by jumping to 
point El here above. 

What happens in the end? If the algorithms returns from all the recursive calls to the 
main algorithm, it is easily seen that energy balanced is reached. Otherwise, we have a 
solution which is decreasing in average energy per sensor (from slice S± towards Sn), and 
which is "locally" energy balanced, for example, we could have: 

±-E x = ]-E 2 = > > Ie 5 = iE 6 > . . . (19) 
Oi h b 3 b A b 5 b 6 

Although not reaching energy balance, we shall prove this solution is optimal. The one 
important thing to observe is the following: If a recursion starts at slice Si, either one of 
the two cases happens: 



The algorithm returns from this recursive call and the solution is locally energy 

h^ 1 h-i 



balanced: —Ei = -^—Ei_i 



• The algorithm does not return from this recursive call and the solution is not energy 
balanced: ^-Ei < £"^7-^-1 • Furthermore, since Si was trying to "catch- up" with Si_i 
and that we set q = 1 (point (fTcjl of the algorithm), it holds that = 0, and thus 
that 

Pi = (20) 
In section |U we use equation (j2TJj) to show that this solution is always optimal. 



3.1.2 Second Case: Little Battery 

The second problem which may occur is when the assumption that all e's are positive, (i.e. 
the assumption from remark |H1 does not hold). From equation (JZJ), we can see that this 
occurs only if some of the slices have little frj's (and thus the title of this subsection). Let 
us first analyse what it means for an e, say e« to be negative. Suppose slice i has, so far, 
ji ejected messages and /j forwarded messages. When it receives k sliding message from 
Si + i, it should eject an fraction to the sink, and pass-on the 1 — e$ rest to the next slice. 
After this, there are ji + kei ejected messages and + k (1 — e,) slided messages. The fact 
that the e's satisfy equation (J2J) ensures that energy balance is conserved (at least locally if 
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we are already into a recursive call as described in section 13. So what is the problem? 
The problem is that since £j is negative, ji becomes negative if 

k ■ ei + ji < <(=> k> ji/ei (21) 

and thus the solution is not physical (a negative amount of messages can not be ejected 
from Si). The fix to this problem has some similarity with the previous one. Whenever a 
slice (say the ith) is about to slide k messages along the network, it should ensure that no 
slice will find itself in a non physical state afterwards by bounding the number of messages 
it allows itself to slide along the network. Suppose a slice Sk wants to slide messages along 
the network, for some fixed k. We call maxSlide the maximum number of messages Sk 
may slide along the network without putting any of its following slices in a non-physical 
state. In order to compute maxSlide, we should remember what happens when k message 
are slided along the network by slice Sf part of it is ejected by each of the slices sliding 
the message, according to the e's, ans therefore only a fc« = k ■ H l j =i+1 (1 — e$) fraction of 
the k initial messages reaches slice Si. maxSlide is defined as the maximum value k such 
that ki€i + ji > for 1 < i < I, or equivalently, the maximum value such that ki < for 
every 1 < i < I such that 6j < 0, and it can be computed by the following procedure. 



Algorithm 3.1: COMPUTEMAxSLlDE(i) 



Input: i, a slice number. 

Output: max, the max number of messages Si can slide. 

global e[ jcomment: This is an array storing the values of the e's 

local F *- 1 



local ejected[ ] 
for (k <— i; k > 1; k 



do | e i ec ^ ec ^[^] *~ F * e[k] 



- F * (1 — e[k]) 
local max <— oo 
for (k i; k > 1; k <— k — 1) 
'local j ejected[k] 
if j<0 
then 

' 3 * j 

if max = oo 
then max <— 



do < 



mi 3 



do < 



else 



return (max) 



{local tmpjnax <— j[k]/j 
if tmpjnax < max 
then max <— tmpjmax 
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Once we have computed maxSlide, we can decide what to do when slice S{ wants to slide k 
messages. If k < maxSlide, then we can simply slide the k messages, but if k > maxSlide, 
we have to be more careful. First, we can partially fulfill the aspiration of Si by allowing 
it to slide maxSlide messages. At this stage, Si still wants to slide k — maxSlide messages 
and one of the previous following slices (say Sk) has a negative e*, and jk = 0. If anymore 
messages are slided, Sk will be in a "non-physical" state. So what we do is that we 
recompute all e^'s for 1 < k < j, 8 but whenever an < and jk = 0, we force e/. to 0. 
What does this do? First of all, it forces the solution to be physical. Second, it breaks the 
relation from equation (j7j) . since for some k's €k is forced to 0. The fact of breaking this 
relation prevents slices from ejecting a negative amount of messages (and thus in some sense 
save some energy), when this would lead them to be in a non-physical state. Thus slice 
Sk will spend more energy than the (locally) energy balanced solution would require, and 
on the other hand, slices following Si (that is S^i to Si) will spend less since the negative 
amount of messages which have been prevented from being ejected where supposed to be 
slided along the network. We are therefore in the presence of a "local peak" , in the sens 
that: 

7~Ek > T E k -\ (22) 

Ok °k-l 

It should be observed that for the rest, energy balance is conserved (at least locally), and 
furthermore whenever such a "peak" appears at Sk, it holds that: 

Pk = 1 (23) 
which is an important fact we shall use to prove that the solution obtained is optimal. 



4 Analysis 

In this section, we prove that our algorithm produces an optimal solution, in the sense 
that it maximizes the lifespan (c.f. definition EJ). 

Convention 11 In this section, we consider a fixed sensor network of size N, with fixed 
event distribution {fl'i}x<i<Ar and fixed battery {bi} 1<i<N . A configuration C of the network 
is the choice of a sliding probability assignment {Pi} 1<i<N for each slice. If C and C are 
two configurations, we use the letters fi and fi to denote the slided messages under the 
configuration C and C respectively. We do the same for the other parameter: ji 's, pi ? s, 
€i ; s and Ei 's. 

Lemma 12 (No Win- Win modification) No configuration is strictly cheaper in terms 
of energy than another configuration: if C and C are two configurations, then there exists 
an i such that: 

j-Ei > —Ei 
bi bi 



8 If we are inside a recursion of the type described in section 13.1.11 we do not recompute all ek 's for 
1 < k < N, but we rather set e s t ar t = 1 and recompute all the e<'s for start < i < N, where start is the 
place where the last recursion took place (c.f. point (|la|) of the enumeration on section 15. 1.11) . 
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Proof. Suppose (absurd) this is not true. Therefore there exists two configurations C and 
C such that 

V* \e, < ^Ei (24) 

and for at least one of the i's 

T E i < T^i (25) 

We now define the following configurations Cq = C and Cn = C, and more generally, Ci 
is the configuration where the i last pis (i.e. pn, Pn-i, ■ ■ -PN-i+i) are the same as the 
Pi's from C, where as the N — i first pi's are the same as the pis from C . Then for each 
1 < i < N, if we use l E and l b to design the energy and battery of configuration Cj, the 
following can be seen to hold: 



^E k =lE k VN>k>N-i (26) 
iT^ (27) 

Where equation (j2EJ) follows from the definition of Ci and where equation (|27jl follows 
from the easy observation that j|r*£?j > ~-Ei combined with equation (|24jl . Next, let 

/c = max < z I < ^Ei \ , which exists by (1251) . Then for every « < the inequality 

I 0i h ) ±<i<N 

in (|2*7j) becomes strict. In particular, for i = (and using the fact that Co = C) it becomes 
^i?o > -^Eq, which is the contradiction we need. □ 
The reason we give this lemma the no win-win modification name is that a principle 
can be derived from it, the no win-win modification principle, which is the following: if 
you have a configuration and you modify it to save energy in some parts of the network, 
then necessarily you will spend more energy in another part of the network. 

In jU], the authors point out that looking at the numerical solutions, one observes that 
an energy-balanced solution mostly uses single-hop data propagation, and only with little 
probability propagates data directly to the sink. The authors then suggest that this is 
an important finding implying that the energy-balanced solution is also energy efficient, 
since it only rarely uses the costly single-hop direct ejection of messages to the sink. Our 
previous lemma enables to easily formalizing this intuition: 



Corollary 13 Any energy-balanced solution is optimal in terms of lifespan: If C is an 
energy-balanced configuration (i.e. Vi = -r^ 1 ), then for every C other configuration, we 

have the following inequality, with equality if and only if C — C : max < max | 
that is, C maximizes the lifespan amongst all possible configurations. 

Next we generalize corollary [T3l 
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Lemma 14 Let C be a configuration of our network. Let max = max j-pj. Let k = 
max < i | p = max \ 9 . Let I = min \i | -P 1 < max \ 10 . The configuration is optimal if and 



bi J 1 bi 

only if the conjunction of the following holds: 

• Pk+i = 

• Pi+i = 1 

Proof. We only give the ideas of the proof. First, notice that slices S& to Si forms a 
tabletop-like maximum of the plotting of slice position against average energy consumption 
per sensor. Since Pk+i = 0, nothing can be done on the left-hand side of the tabletop to 
lower it. Second, since pi + i = 1, the tabletop can not rely on the slices on its right to 
take-on a larger part of the message sliding towards the sink. The only solution to produce 
a better solution than C (i.e. if C was not an optimal solution), would therefore be to 
modify the probabilities from pk to pi+i, i.e. to reorganize the configuration "inside the 
tabletop". The final point is to notice that this will break the energy balance of the 
tabletop, increasing the maximum (using the principle of no win-win modifications from 
lemma IP2|). □ 



Theorem 15 Our algorithm always produces an optimal solution 

Proof. We would be done if we could prove that our algorithm always produces a solution 
where the maximal is reached at a tabletop with pi = on the left and pi = 1 on the right, 
since this enables the call to lemma ITU To see that this is the case, the main ingredients 
are equations (ffilj) and (f2l3|) . We leave the easy details to the reader. □ 
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A Pseudocode of the Algorithm 

We provide hereunder a pseudode of the algorithm presented in this paper. This pseu- 
docode is rigourously based on a perl implementation of the algorithm which was validated 
on various inputs. 
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Algorithm A.l: COMPUTEOptimal(A^; g[ ]; b[ ]; d[ ]) 



global /[ ]; j[ ]; max = oo; recLevel = 0; startPosition = 1; e[ ] 
comment: Following arrays used as stacks while changing recLevels 

global start Positions[ };maxs[ ]; epsi/ons[ ] 
main 

global i = 0; initialG[ ] = <?[]; 
for i <- 1 to JV 

do {/[i] = = 
e[ ] <— EPSILONS(l) 
while i < N 

<i <— i + 1 
PUSH(recLeve/s[ ], recLevel) 
if i = 1 

comment: First step 
then 

<- 



do 



else 

comment: From second step to the Nth 

f local El <- /[i - 1] + j[i - 1] * d[i - l] 2 
local idealj <— avgNrj(? - 1) * o[i]/d[i] 2 
if idealj > g[i] 

comment: Not enough messages to stay at this level 
then 

do f EJECT((?[i]) 

\downOneLevel() 



else 

comment: Enough messages to stay at this level 

' EJECT(idealj) 
while (g[i] > 0) 

'local nrj Delta <— max — AVGNRj(i) 
local msgToGoUp <- t . rr^fj^ r . m 

if recLevel = or g[i] < msgToGoUp 
comment: Slide the rest 



do < 



do < 



do < 



then SLlDE(#[i]) 



else 



comment: Slide enough to go up one level 
J siMsimsgToGoUp) 



do 



1 upOneLevelQ 



procedure AVGNRj(pos) 
return (iH)(/M+j'H * d[pos} 2 )) 



procedure EPSILONS(/irst; option) 
local e[ ] 

for k <— 1 to /irsi 

do {e[ib] = 1 
for k <— first + 1 to N 
'local A 

b[k] 

local 5 <- (e[ife - 1] * (d[k - l] 2 - 1) + l)/b[k - 1] 
e[k] «- 

; if option = "withCaution" 
then 

{if e[Jfe] < and j[Jfe] = and fc < % 
then 
do {e[ib] = 
return (epsilons[ ]) 



procedure upOneLevel() 
local max <— POP(maa;s[ ]) 
startPosition <— POP(startPositions[ ]) 
local imp[ ] <— POP (epsi/ons[ ]) 
e[ ] <— tmp[ ] 

recLevel <— recLevel — 1 



procedure downOneLevel() 
P\JSK(maxs[ },max) comment: store old max 

max <— avgNrj(« - 1) 
P\JSu(startPositions[ ], startPosition) 
startPosition <— i 
local tmp Array [ ] = e[ ] 

PUSH(epsi/ons[ ],tmpArray[ ]) comment: store old epsilons 

e[ ] <— EPSILONS(i) 
recLevel <— recLevel + 1 



procedure EJECT(eject) 
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j[i] <- eject 
g[i\ <- 

procedure SLIDE(F) 
slideCarefull(F) 

procedure slideCareless(F) 

<7[i] <- <7[i] - F 

for (A; <— A; > 1; k <— A; — 1) 

f /[*]<-/[*]+ ^ * (1 " e[*]) 
do < j[jfc] <- j[k] +F*e[k] 

(F <- F * (1 - e[k]) 

procedure slideCarefull(F) 
local maxCarelessSlide <— computeMaxSlide() 
if maxCarelessSlide = oo or F <= maxCarelessSlide 
then 

do (slideCareless(F) 
else 

{SLlDECARELESS(maxCareless Slide) 
F <— F — maxCarelessSlide 
e\ ] <— EPSILO~NS(startPosition,"withCaution") 
SLIDE(F) 

procedure computeMaxSlideQ 
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F <— 1 comment: Simulated slide of one packet from the current pos ( i ) 

local ejected[ ] 
for (k <— i; k > 1; k <— k — 1) 
J ejected [A;] <— F * e[k] 
d ° \F <- F(l - e[k\) 
local max <— oo 
for (k <— i;k > 1; <— — 1) 
'local j <— ejectedfA;] 
if J < 
then 

j < — j; 

if (max = oo) 
then 

do max <— ^ 
else 

{local tmp m ax <— ^jp 
if tmp m ax < max 
then 
do max = tmp m ax 
return (max) 



do < 



do < 
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